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The coupled equations that describe the effect of large-scale magnetic and velocity fields on forced high-diffusivity 
magnetohydrodynamic flows are investigated through an extension of mean field electrodynamics. Our results generalise 
those of Radler & Brandenburg (2010), who consider a similar situation but assume that the effect of the Lorentz force on 
the momentum equation can be neglected. New mean coupling terms are shown to appear, which can lead to large-scale 
growth of magnetic and velocity fields even when the usual Q-effects are absent. 
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1 Introduction 

Magnetohydrodynamic turbulence, particularly in the con- 
text of astrophysics, is often studied within the framework 
of mean field electrodynamics, an elegant closure theory 
that describes the evolution of mean quantities in terms of 
transport coefficients determined from averaged small-scale 
properties of the flow and magnetic field (see, for example, 
Moffatt 1978, Krause & Radler 1980). 

The simplest form of the theory considers the case when 
the magnetic field B can be regarded as kinematic, evolv- 
ing under the influence of a velocity U, but exerting no in- 
fluence back on the velocity. Only the induction equation is 
then of significance. The theory proceeds by decomposing 
the velocity and the magnetic field into mean and fluctuating 
parts, 

U = U + u, B = B + b, (1) 

under an averaging procedure obeying the Reynolds rules. 
The most common practice is to adopt spatial averaging, as- 
suming a distinct scale separation between that of the large- 
scale (mean) fields and the small-scale (fluctuating) fields. 
If, for simplicity, we assume that there is no mean flow (i.e. 
Uq = 0), then averaging the induction equation leads to the 
following equations for the mean and fluctuating magnetic 
fields, 

dB 



dt 



= V x S + tjV'Bq, 



db 

di 



V x (u x B ) + V x G + i]V 2 b, 



(2) 



(3) 



where £ = (u x b) is the mean electromotive force (emf), 
G = (u x b) — (u x 6) and r\ is the magnetic diffusivity. 
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Equation (0 can be expressed as 

Lib) = V x {u x B ) , 



(4) 



where £ is a linear operator. If the small-scale field can grow 
in the absence of a mean field, i.e. if the equation C(b) = 
has exponentially growing solutions, then a small-scale dy- 
namo is operative, and the interpretation of any mean field 
equations becomes problematic (Cattaneo & Hughes 2009). 
If, however, non-decaying solutions for b depend crucially 
on a non-zero mean field Bo, then b, and hence £, depend 
linearly on B, and it is customary to express £ as a series 
in B and its spatial derivatives (though see Hughes & Proc- 
tor (2010) for a discussion of the significance of the omitted 
temporal derivatives), usually written as 



a>ijB j 



a 9B 0j 
' Pijk-^ r 

axk 



(5) 



where it is anticipated that the large spatial scale of Bq will 
lead to rapid convergence. The tensors ay and /3y fe depend 
on the statistical properties of the velocity field and on rj. 

The above kinematic formulation can be regarded as 
the study of the evolution of weak, long wavelength mag- 
netic field perturbations to a pre-existing non-magnetic flow. 
One may then extend this idea so as to consider long wave- 
length perturbations, in both the magnetic field and the ve- 
locity, of a magnetohydrodynamic [MHD] state involving 
both small-scale field and flow; here the perturbations are 
again infinitesimal, but the background state is, in general, 
fully nonlinear. This problem has been tackled recently in 
a variety of ways. Courvoisier, Hughes & Proctor (2010a) 
considered the evolution of linear, three-dimensional per- 
turbations to two-dimensional MHD basic states resulting 
from prescribed forcings in the momentum equation in the 
presence of a background magnetic field with non-zero flux. 
Courvoisier et al. (2010a) stressed the importance of treat- 
ing the magnetic and velocity fields on an equal footing, and 
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showed how the linear evolution of the mean field and mean 
flow are both dependent on terms proportional to both the 
mean field and the mean flow. In an Appendix, Courvoisier 
et al. (2010a) outlined the formal extension of the theory to 
describe perturbations of a three-dimensional, MHD back- 
ground state. This work is extended in Courvoisier et al. 
(2010b), which considers in more detail the evolution of 
long wavelength perturbations to a fully nonlinear three- 
dimensional, small-scale turbulent state, as may result from 
the nonlinear saturation of a small-scale dynamo. Radler & 
Brandenburg (2010) have also recently considered the na- 
ture of the emf in MHD turbulence. They made the simpli- 
fying assumption that both the basic state and perturbations 
to it were kinematic, omitting the Lorentz force throughout. 
However, they allowed a basic MHD state to arise by pre- 
scribing an external electromotive force in the induction eq- 
uation, in addition to an external body force in the momen- 
tum equation. They then considered the influence of a weak 
mean velocity (but no mean magnetic field), their main re- 
sult being that the ensuing mean emf has a component pro- 
portional to the mean flow. This result had also been noted 
by Courvoisier et al. (2010a: equations (A. 17) and (A. 18)). 

It is interesting to note that in the model studied by 
Radler & Brandenburg (2010), it is not actually necessary to 
make the simplifying assumption of neglecting the Lorentz 
force. In this paper we therefore extend their analysis by 
retaining the Lorentz force throughout, thus allowing us to 
study the evolution of more general MHD states. However, 
we ignore rotation and any effects of mean field gradients, in 
order to focus on the most basic effects: the mean emfs and 
Reynolds and Maxwell stresses that arise owing to imposed 
mean velocity and magnetic fields. We adopt a slightly dif- 
ferent approach to that of Radler & Brandenburg (2010), 
making the assumption that the magnetic Reynolds number 
Rm is small, so as to make analytic progress. In order that 
we can consider a basic MHD state for small Rm, we re- 
tain the external emf in the induction equation, since, with- 
out this, any small-scale field would simply decay. We show 
that when considering a background MHD state (as opposed 
to a purely hydrodynamic state), the magnetic and velocity 
perturbations must, for consistency, be treated on the same 
footing. Thus, one must consider not just the mean emf in 
the averaged induction equation, but also the mean stress 
tensors in the averaged momentum equation. 

2 The mean emf and the mean stress tensor 

Our treatment starts with a basic MHD state consisting of 
fluctuating, small-scale velocity and magnetic fields, U and 
B, that are driven by an applied body force and an applied 
electromagnetic forcing. The fluid has kinematic viscosity 
v and magnetic diffusivity 77. The basic state is described by 
the following non-dimensional equations, 

R m { 9 -£ + U.Vu) (6) 

R m (-vn + X B ■ VB) + P m V 2 U + f, 



dB 

+ U ■ MB I 
= Rm (B ■ VU) + V 2 J3 + h, 



/,',„ 1 — + u ■ vsj 



(7) 



where n represents the total pressure. The velocity is scaled 
with U = J-C 2 jr\ and the magnetic field with B = R.C? jr\, 
where T, H and C are representative scales for the applied 
body force, the applied electromagnetic force and length. 
The dimensionless parameters are the magnetic Reynolds 
number R m = UL/rj = J-£ 3 /rj 2 ; the magnetic Prandtl 
number, P m = v/r}\ and \ = (H 2 1 PHo) /J 72 , where p is the 
(constant) density and /jq the magnetic permeability. The 
parameter x, which may be written as M.^ 2 , where M.a 
is the Alfvenic Mach number, describes the relative impor- 
tance of the magnetic field and the velocity in the basic state. 
The situation investigated by Radler & Brandenburg (2010) 
is recovered by letting \ — > 0. 

We now introduce applied uniform velocity and mag- 
netic fields, Uq and Bo, with the aim of evaluating the 
essential mean quantities in the momentum and induction 
equations. We assume that the energies in the mean flow 
and field are of the same order, or smaller, than those of the 
basic state. If we express the total velocity as u+Uq and the 
total magnetic field as b + Bo, then the governing equations 
take the form, 



Rn 



du 
~dt 



u ■ Vu + Uq ■ Vu 



R rn (-Vtt + X b ■ Vb + x-Bo • V6) 
+ P m V 2 u + f, 



, db 

/?,„ ( — + u-Vb + U -Vb 



R m (b ■ Vu + B Q ■ Vu) + V 2 b + h. 



(8) 



(9) 



We suppose that both R m and the Reynolds number R e = 
PmRm are small, with P m of order unity. We may thus ex- 
pand u and b in powers of R m , namely 

u = u<°> + R m u^ + ■■■ , b = &<°> + R m b^ + ■■■ . 

(10) 

Similarly, we may expand the fields u' = u — U and 
b 1 = b B, induced by the imposed mean velocity and 
magnetic fields, as 

u' = R m u' {1) + ■■■ , b' = R m b' (1) + ■ • • , (11) 

where 

P m VV« - (U ■ W°) - X Bo ■ V&<°>), (12) 

V 2 b'« = (U ■ V&(°> - B • (13) 

The right hand side of ( fT2l is solenoidal, and hence there is 
no pressure gradient term. We assume that there is no mean 
emf in the basic state, described by equations (O and (0; 
hence the leading order contribution to the emf is given by 



(14) 
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To aid calculation, we now further suppose that / and h are 
independent of t, spatially periodic and monochromatic, so 
that 

V 2 / = -A 2 /, W 2 h=-\ 2 h, (15) 

with |A| ~ Cr 1 . It can be shown that it' ) and are then 
also monochromatic. Hence we can write 

(u<°> x = -A- 2 (VV°> x 6'«) (16) 
= -A-V°> x V 2 6' (1) ), (17) 

and similarly for the other term in ( TBI . All surface terms in- 
volved in the integrations by parts needed to transform ( fTol i 
to ( TP71 ) vanish, by assumption of spatial periodicity. Then 
using equations Sl2\ and ( TT~3T > we obtain exactly, 

A 2 i? m ^ = u ■ {(l - p-^vufh^) 

+ B,.{e ijk {ufWuf + P -\Wbfbf)). 

(18) 

The terms in Bq may be recognised as the usual kinetic and 
magnetic a-effects, while the term in [To has the same form 
as that discussed by Radler & Brandenburg (2010) in a more 
general context. 

Thus far there is little that is new. However, crucially, we 
can see that because the momentum equation now contains 
Bo, through the Lorentz force, we also have a non-trivial 
expression for the mean stress tensor, R^, at leading order, 
namely 

iC^H< (1 M 0) +^; (1) > 

Under the assumptions ( [TBI ), and via manipulations analo- 
gous to those leading to ( fT8l , we obtain 

X 2 R^R r] = xB ■ ((1 + P^XV&f uf + Vbfuf)). 

(20) 

It might have been expected that there would also be a term 
proportional to Uq in this expression. Were such a term to 
exist, it would be analogous to the AKA (anisotropic ki- 
netic a) effect introduced by Frisch, She & Sulem (1987). 
There are such contributions present here, but they cancel 
out; in fact it can be shown, in general and not just for 
monochromatic velocity fields, that the AKA effect van- 
ishes in the first order smoothing approximation, for steady 
forcing as assumed. However, a non-zero contribution may 
be expected for more general forcings. In spite of this, it is 
plain from the form of the terms in expressions ( TT~8T > and ( f20b 
that the evolution equations for the large-scale velocity and 
magnetic fields will generic ally be coupled together. Inter- 
estingly, both types of coupling term depend on the tensor 

<2I) 

The above calculations have been performed for uni- 
form, steady Uo and B . However, and this is the under- 
lying principle of mean field electrodynamics, the results 



for the mean emf and the mean stress tensor can be used to 
determine the evolution of velocity and magnetic fields de- 
pending slowly on space and time. Even for small Rm, if the 
length scales are sufficiently long then the terms involving 
first derivatives dominate over diffusion terms, which may 
therefore be neglected. In this case, the mean momentum 
and induction equations take the form 

d T U 0l + d X] Rij = -d Xt P , (22) 

drBoi = eijkdxj£k, (23) 

where T and X are long time and space variables. We now 
ignore the familiar a-effect terms, and look for solutions 
proportional to e zK - x + sT . After appropriate scaling of the 
time and space variables, to absorb the factors of X 2 R^ n 1 , 
and elimination of the pressure using V • Uq — 0, we obtain 
(writing Wf jk = Q ikj ± Q klj ) 

sB 0i - (1 - P^) W M ■ iK k U ol , (24) 
sUoi - (1 + P" 1 ) xSijWf kl ■ iK k B ol , (25) 

where Sij = 5y — KiKj/\K\ 2 . While W + , being symmet- 
ric in its first two arguments, will vanish in isotropic situa- 
tions, there is no reason for it to vanish in general. Indeed, 
for a simplified problem involving a 2D basic state, Cour- 
voisier et al. (2010a) have exhibited cases where this term 
is non-zero. Because the equations are coupled, exponential 
growth is possible even without the a-effect. The square of 
the growth rate s is given as the eigenvalue of the matrix 

- X (l - Wr kl K k K p S lm W+ pg , (26) 

and it appears that solutions with a positive real part of s 
can be found in a wide variety of situations, provided that 

Pm + 1. 

3 Discussion 

Starting from a small-scale, MHD basic state, forced via 
a body force and an externally applied emf, we have con- 
sidered the mean emf, £, and the mean Reynolds-Maxwell 
stress tensor, Rij, resulting from the imposition of spatially 
uniform, steady, velocity and magnetic fields. Under the as- 
sumption that the fluids and magnetic Reynolds numbers 
are small, we have obtained analytic expressions for £ and 
Rij, ( fT8l and ( f20b . These are linear in the mean fields, with 
transport coefficients dependent on the means of various 
quadratic fluctuating quantities of the basic state. 

The main point of our paper is that in a true MHD state, 
in which neither the velocity nor the magnetic field can 
be regarded as dominant, it is imperative that the velocity 
and magnetic fields are treated on an equal footing. This 
is in sharp contrast to traditional mean field electrodynam- 
ics, which is essentially a kinematic theory, describing the 
evolution of a (formally weak) large-scale magnetic field 
under the influence of a prescribed velocity. It is then only 
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the mean emf that is of interest; furthermore, the mean emf 
depends only on the fluctuating velocity field. In the case 
considered in this paper, however, £ and Rij are of equal 
importance and, moreover, they both depend on both the 
fluctuating velocity and the fluctuating magnetic field. We 
believe that this is an important message, in that it highlights 
the shortcomings in, for example, considering only the emf, 
and then only the part proportional to the mean magnetic 
field (the a-effect), for MHD turbulence. Although it is of- 
ten claimed that the a-effect in the nonlinear (MHD) regime 
can be expressed as the difference between the flow helicity 
and the current helicity (essentially the Bq term in (fTSI)). our 
analysis shows clearly that this tells only part of the story. 
In the general case, the coupling between the mean velocity 
and magnetic fields is such that the growth of a large-scale 
disturbance can be expected even when the usual kinematic 
and magnetic a-effects vanish. Such a mode of instability 
is not accessible if the influence of the Lorentz force on the 
fluid momentum is dropped. 
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